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We study the statistics and characteristics of rare intense events in two types of two dimensional 
Complex Ginzburg-Landau (CGL) equation based models. Our numerical simulations show finite 
amplitude collapse-like solutions which approach the infinite amplitude solutions of the nonlinear 
Schrodinger (NLS) equation in an appropriate parameter regime. We also determine the probability 
distribution function (PDF) of the amplitude of the CGL solutions, which is found to be approx- 
imately described by a stretched exponential distribution, P(jA|) ~ e~' A , where r\ < 1. This 
non-Gaussian PDF is explained by the nonlinear characteristics of individual bursts combined with 
the statistics of bursts. Our results suggest a general picture in which an incoherent background 
of weakly interacting waves, occasionally, 'by chance', initiates intense, coherent, self-reinforcing, 
highly nonlinear events. 

PACS numbers: 02.30.Jr, 03.65.Ge, 05.45.-a, 52.35.Mw 
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I. INTRODUCTION 

Many spatio-temporal dynamical systems show rare in- 
tense events. One example is that of large height rogue 
ocean waves [jp. Another example occurs in recent ex- 
periments on parametrically forced surface waves on wa- 
ter in which high spikes (bursts) on the free surface form 
intermittently in space and time jl . Other diverse physi- 
cal examples also exist (e.g., tornados, large earthquakes, 
etc.). The characteristic feature of rare intense events is 
an enhanced tail in the event size probability distribution 
function. Here, by enhanced we mean that the event size 
probability distribution function approached zero with 
increasing event size much more slowly than is the case 
for a Gaussian distribution. Thus these events, although 
rare, can be much more common than an expectation 
based on Gaussian statistics would indicate. The central 
limit theorem implies Gaussian behavior for a quantity 
that results from the linear superposition of many ran- 
dom independent contributions. Non-Gaussian tail be- 
havior can result from strong nonlinearity of the events, 
and enhancement of the event size tail might be expected 
if large amplitudes are nonlinearly self-reinforcing. Such 
nonlinear self-reinforcements is present in the nonlinear 
Schrodinger (NLS) equation. In particular, the two di- 
mensional NLS equation, 



dA 
~dt 



-ia\A\ 2 A - if3 \/ 2 A. 



(1) 



exhibits "collapse" when the coefficients a and (3 have the 
same sign ||. In a collapsing NLS solution the complex 
field approaches infinity at some point in space, and this 
singularity occurs at finite time. The NLS is conservative 
in that it can be derived from a Hamiltonian, dA/dt = 
-iSH/SA*, where H[A, A*] = ± f[a\A\ 4 + /3| y A\ 2 }dx. 
In the case of nonconservative dynamics, inclusion of 



lowest order dissipation and instability terms leads to 
the complex Ginzburg-Landau (CGL) equation [Q. The 
CGL equation has been studied as a model for such di- 
verse situation as chemical reaction ||, Poiseuille flow 
H , Rayleigh-Bernard convection Q , and Taylor-Couette 
flow ||. In the limit of zero dissipation/instability the 
CGL equation approaches the NLS equation. For small 
nonzero dissipation/instability, the CGL equation dis- 
plays a solution similar to the NLS collapse solution, but 
with a large finite (rather than infinite) amplitude at the 
collapse time ||. Furthermore, over a sufficiently large 
spatial domain, these events occur intermittently in space 
and time. Thus, in this parameter regime, the CGL equa- 
tion may be considered as a model for the occurrence of 
rare intense events. 

In this paper we study the statistics and characteristics 
of rare intense events in a two-dimensional CGL-based 
model. The probability distribution function (PDF) 
of the amplitude of the solutions is observed to be 
non-Gaussian in our numerical experiments. This non- 
Gaussian PDF is explained by the nonlinear characteris- 
tics of individual bursts combined with the statistics of 
bursts. The model equation we investigate is 



r) A 

— = ±A-(l+ia)\A\ 2 A+(l-i/3) V 2 A 
+(S r + i8i)A*, 



(2) 



where (5 r + iSi)A* is a parametric forcing term |J. We 
will consider two cases: one without parametric forcing 
(S r = 5i = 0) in which case the plus sign is chosen in 
front of the first term on the right-hand side of (^) (Eq. 
(@) is then the usual CGL equation), and one with para- 
metric forcing, in which case the minus sign is chosen. As 
previously discussed, we choose our parameters so that 
our model, Eq. (|J), is formally close to the NLS equation 
([0). That is, we take a, (3 S> 1, S r , §i, and for our numeri- 
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FIG. 1: Solutions of the CGL model, (a) Snapshot of the 
amplitude \A\ for L = 20tt, a = f3 = 30, At = 10~ 5 ,6 r = 
Si — and a 256 x 256 grid, (b) Amplitude profile versus 
time. Solid line indicates |j4.| ma:c of the whole system, while 
thick solid lines indicate maximum amplitude of the local- 
ized events(" bursts"). The dashed line indicates the average 
amplitude of \A\ over the whole system |A[ a „ 9 ~ 0.3. 



cal solutions we will restrict attention to the case a = (3. 
Note that the coefficient ±1 for the first term, as well 
as the ones in (1 + id) and (1 + i(3) represent no loss of 
generality, since these can be obtained by suitable scaling 
of the time(i), the dependent variable(A), and the spa- 
tial variable(x). In Section II, we discuss the amplitude 
statistics of our two-dimensional CGL models with and 
without the parametric forcing term. We find that the 
PDFs are approximately described by a stretched expo- 
nential distribution, P(|A|) ~ exp(— \A\ V ), where r\ is less 
than 1. In Section III, we investigate the characteristics 
of individual bursts. We compare our numerical CGL 
results with known collapse solutions of the NLS equa- 
tion. The maximum amplitude obtained by many bursts 
(or the "event size" statistics) is discussed in Section IV. 
Section V presents further discussion and conclusions. 

Our results lead us to the following picture for the 
occurrence of rare intense events in our system. Linear 
instability and nonlinear wave saturation lead to an inco- 
herent background of small amplitude waves. This back- 
ground is responsible for the observed small \A\ Gaus- 
sian behavior of our probability distribution functions. 
When, by chance, the weakly interacting waves locally 
superpose to create conditions enabling nonlinear, coher- 
ent self-reinforce, a localized, collapse-like event is initi- 
ated. Collapse takes over, promoting large, rapid growth 
and spiking of A. This is followed by a burn-out phase 
in which the energy in rapidly dissipated due to the gen- 
eration of small scale structure by the spike. We believe 
that elements of the above general picture may be rele- 
vant to a variety of physical situations where rare intense 
events occur (e.g., the parametrically driven water wave 
experiments in Ref. ^]). 



II. AMPLITUDE STATISTICS 
A. 2D model without a parametric force 

(S r =5i=0) 

We first consider Eq. <^) with S r = Si = and with the 
plus sign in the first term on the right hand side of the 
equation. Figure |l](a) shows a representation of |A(x, t)\ 
[from numerical solution of Eq. (g)] at a fixed instant t, 
where large values are indicated by darker grey shades. 
As a function of time, the localized dark shades occur 
in an seemingly random manner, become darker (i.e., in- 
crease their amplitude) and then go away (become light). 
Furthermore, the maximum amplitudes also display ap- 
parent randomness, [see Fig. [j](b)]. As shown in the next 
section (Sec. Ill), although the occurrence of these in- 
tense events is apparently erratic in time and space, indi- 
vidually these events are highly coherent. In this section, 
we will study the statistical properties of A(x, t). 

Our numerical solutions of (Q) employ periodic bound- 
ary conditions on a 256 x 256 grid. We choose the pa- 
rameters, a and (3, large enough (a = (3 = 30) so that 
the solutions of our model are close to solutions of the 
NLS equation. We choose the time step small enough 
to satisfy the condition for unconditional stability of our 
second-order accurate time integration(At = 10~ 5 ). We 
use random initial condition (at t = 0, we generate ran- 
dom values for amplitudes and phases at each grid point). 
Localized structures, "bursts", develop very rapidly and 
occur throughout the spatial domain. The typical life 
time of a burst (St) is approximately 0.2 time units. The 
maximum amplitudes of bursts are different for different 
burst events. 

Imagining that we choose a space-time point (x, t) at 
random, we now consider the probability distribution 
functions for \A\ (the magnitude of A), A r = Re[A] 
(the real part of A), and Ai = Im [A] (the imaginary 
part of A). We denote these distribution functions 
P(\ A\), P r (A r ), Pi(Ai), and we compute them via his- 
togram approximations using the values of |A|,A r , and 
Ai from each of the 256 x 256 grid points at many time 
frames |l0| . We find that these distributions are indepen- 
dent of the periodicity length L used in the computation 
as long as it is sufficiently large compared to the spatial 
size of a burst, but is not so large that spatial resolu- 
tion on our 256 x 256 grid becomes problematic. In our 
computations of P,P r and Pi, we choose L = 207r. 

Figures 2 show the numerically computed probability 
distributions P r (A r ) [Fig. 2(a)] and P l (A l ) [Fig. 2(b)] 
plotted as open circles. Since Eq. (||) with S r = Si = is 
invariant to the transformation A — > Aexp(i(j)) (where <f> 
is an arbitrary constant), we expect the distribution P r 
and Pi to be the same to within the statistical accuracy 
of their determinations. This is born out by Figs. 2. 
In order to highlight the essential role that nonlinearity 
plays in determining these distribution functions, we have 
also computed them after randomizing the phases of each 
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FIG. 2: Probability distribution functions obtained from nu- 
merical solution of Eq.(2) using the same parameters as in 
Fig. 1. The circles are data for P r and Pi, while the pluses 
are the probability distributions P' r and P[ obtained from the 
phase randomized amplitude. 




FIG. 3: Probability distribution functions before randomiz- 
ing the phases of the solutions (o) and after randomizing the 
phases (+). Note the horizontal axis is \A\* , where the expo- 
nent rj = 0.8 is chosen to yield approximately linear depen- 
dence of log[P(|A|)] on \A\ V for large values of \A\. 



Fourier component. That is, representing ^4(x, t) as 

A(x,t) = ^a k (t)exp(zk-r), (3) 

k 

where k = (2mir/L, 2mr/L), we form a new amplitude, 
A'(x,t) as 

A'(x,t) = ^2a k {t)exp(ik-r + i8 k ), (4) 

k 

where for each k, the angle 8^ is chosen randomly with 
uniform probability density in [0,27r]. The probability 
distribution functions for the real and imaginary parts 
of the randomized amplitudes A' are shown as pluses 
in Figs. H Note that by construction, A and A' have 
identical wavenumber power spectra. Due to the random 
phases, A' at any given point x can be viewed as a sum of 
many independent random numbers (the Fourier compo- 
nents). Hence the P r and Pi distributions are expected 
to be Gaussian, logP r! ,; ~ [(const.) — (const.) A^ i ]. This 
is confirmed by the semi-log plots of Figs. 2, where the 
data plotted as pluses can be well-fit by parabolae. 

The above comparisons with the phase randomized 
variable A' are motivated by imagining the hypothetical 
situation where the amplitude is formed by the superpo- 
sition of many noninteracting linear plane waves. In that 
case we would have an amplitude field of form 

^ bk(t) exp(zk ■ r + iuj^t). (5) 

k 

Because Wk is different for different k, the phases become 
uncorrelated for sufficiently large time t. 

Comparing the data from A and A' in Figs. 2, for 
small values of A r and Ai, the PDFs are nearly Gaus- 
sian. This can be attributed to near linear behavior of 
small amplitude waves. On the other hand, for the tails of 
the distributions, we note substantial enhancement rela- 
tive to the Gaussian distributions. These must be due to 
coherent phase correlations resulting from nonlinear in- 
teraction of different wavenumber components of A. Such 



phase coherence is implied by the observed coherent lo- 
calized burst structures. 

Figure 3 shows the numerically obtained distribution 
P(|A|) plotted as circles and the probability distribu- 
tion for the phase randomized amplitude \A'\ plotted as 
pluses. Again, the enhancement of the large amplitude 
tail is evident. Note that the vertical axis in Fig. 3 is 
logarithmic, while the horizontal axis is \A\ n . Here we 
choose the power n = 0.8 so that the large \A\ data in 
this plot are most nearly fit by a straight line. That is, 
we attempt to fit P(|^4|) using a stretched exponential. 
The slope of the dashed straight line in the figure is cho- 
sen to match the large |^4| data. Thus, over the range of 
\A\ accessible to over numerical experiment, we find that 
the enhanced large |^4| tail probability density is roughly 
fit by a stretched exponential, 

P(L4|) ~ exp(-C|A| - 8 ). (6) 

B. 2D model with parametric forcing (5 r ,5i ^ 0) 

We now report similar results for the case of paramet- 
ric forcing, Eq. (|^) with S r ,5, ^ and the minus sign 
chosen in the first term on the right side of @. In this 
case, instability of small amplitude waves is caused by the 
parametric forcing (nonzero 8 r ^) and the —A term rep- 
resents a linear wavenumber independent damping. This 
model for parametric forcing was introduced pd| ] and has 
been used to model various situations. One such situa- 
tion is that of periodically forced chemical reactions [jl2| . 
Our motivation for considering this model is the Fara- 
day experiments on strong parametric forcing of surface 
water waves in Ref. In that work intermittent for- 
mation of large localized surface perturbations results in 
splash and droplet formation. 

Parameters for our numerical simulations are the same 
as in Sec. IIA except that now S r = Si = 5. Again 
coherent localized structures, "bursts", develop rapidly 
and occur intermittently throughout the spatial domain, 
Fig. |^(a). As in Sec. IIA, the typical life time of a burst 
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FIG. 4: Solutions of the model with parametric forcing, (a) 
snapshot of |j4|. Dark regions have high amplitudes, (b) 
P(|A|) versus \A\ V , where r\ = 0.8. (See captions in Fig. [|) 

is less than 0.2 time units, and the maximum amplitudes 
of bursts are different for different bursts. 

The PDF, P(|A|) again shows a stretched exponential 
tail with exponent rj = 0.8, Fig. ^(b). The circles in Figs. 
H show the PDFs of the real and imaginary parts of A, 
while the pluses are data for the PDFs after randomiz- 
ing the phase. The shape of the PDFs around the central 
part are nearly Gaussian. In contrast, at large amplitude 
the PDFs are significantly non-Gaussian. A major differ- 
ence with the case S r = Si = is that, with parametric 
forcing, the model is not invariant to A — > Ae 1 ^, and thus 
P r and Pi may be expected to evidence differences not 
present for 6 r = Si = 0. This is seen to be the case in 
Figs. H 



III. CHARACTERISTICS OF INDIVIDUAL 
BURST EVENTS 

Solutions of the CGL equation with large a and (3 may 
be expected to have feature's in common with solutions 
of the NLS equation. It is known that the NLS equation 
yields localized events which develop finite time singular- 
ities where the amplitude becomes infinite B . While it 




FIG. 5: P r {A r ) and Pi(Ai) with parametric forcing. The 
numerical parameters are the same as those for Fig. [4| The 
circles are PDFs for A r and Ai before randomizing the phases 
of the solutions, while pluses are PDFs after randomizing the 
phases. 
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FIG. 6: Self-similar bursts, (a) Enlarged plots of a burst 
at ti(O) = 10.448, at grid point (x,y) = (15.2,22.1), and 
i 2 (A) = 10.530, at grid point (51.1, 25.5), and t 3 (D) = 10.644 
at grid point (57.4,38.8). (b) Scaled profiles, where |j4| = 
|-A|/L, x = x/L, and L — \A\ m ax at ti, t2. and t%. The solid 
line represents the radial solution of Eq. (H) . 

is difficult to understand the dynamics of the solutions 
of CGL equation from direct rigorous analysis, the solu- 
tions of the NLS equation are relatively well understood. 
Thus, we analyze the dynamics of individual CGL bursts 
guided by the known localized self-similar collapsing so- 
lution of the NLS equation. 

The NLS equation has a special solution (p| of the 
form 

A = e m R(r), r = v^ + y 2 , (7) 
where the radial function R(r) satisfies 

(£ + i|)*-t« + rf-o, m 

^ =0, fl(oo)=0, 

where a = j3, £ = 9//3 . Since (|l|) is invariant under the 
scaling transformation 113], 

A(x, t) — » L^)- 1 A(k, T)e lii|K|2/4 , (9) 
where Lit) tends to zero as t* — ► t, t < t* and 

K = wr T = Lim ds ' (10) 

a family of collapsing solutions of the NLS is given by 
the rescaled solution of Eq. (||) . 

With these considerations, we test the expected ap- 
proximate self-similarity of individual bursts observed in 
our numerical solutions of Eq. We consider three 

typical bursts that occur at different times and spatial 
positions. In particular, we choose these three as the 
three dark regions in Fig. |l|(a) whose spatial maxima 
have the time dependences shown as thick solid lines in 
Fig. 0(b). 

In Fig. ||(a) we plot the ^-dependence of |A| at con- 
stant y for each of these bursts at the time that they 
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FIG. 7: Self-similarity of single burst, (a) Enlarged plots of 
a burst at ti(Q) = 10.50, t 2 (D) = 10.52, t 3 (A) = 10.53, 
<4(o) = 10.54, and ts(x) = 10.55. (b) Scaled profiles, where 
[j4| = \A\/L, x = x/L, and L = \A(r)\ ma x at tp t 2 , t s , and tn. 
The solid line represents the solution of Eq. (| 



reach their maximum amplitude (the positions in x of 
the maxima have been shifted to x = and the con- 
stant y value for each is at the location of |v4.| moa; ). Note 
that, when they reach their maxima, the three bursts 
have different amplitudes and width. We rescale the am- 
plitude and spatial coordinate as suggested by (0) and 
(|l0|), \A\ = \A\/L and x = x/L, and we take L = \A\ max 
(which normalizes |A| maa; to one). The resulting data are 
plotted in Fig. |6|(b) along with the solution of Eq. (||). 
[We again rescale R(r) using (g) and @, and we note 
that, after this rescaling, the result is independent of the 
value of £ = 0.1 in Eq. (0).] The three burst profiles 
show evidence of collapsing onto the theoretical radial 
solution. 

Now, we consider the time dependence of a single 
burst. We select the burst at the grid point (x,y) — 
(51.1,25.5) (see caption to Fig. 0) and investigate the 
evolution of its shape and height. We display profiles of 
the burst at 5 different times in Fig. 0(a). Rescaling each 
profile using Eq. (0) and Eq. ([To]) , and defining L in the 
same way as before, the four profiles at the first four times 
approximately collapse onto the radial solution of (0) as 
shown in Fig. 0(b). When the burst reaches its max- 
imum amplitude, the amplitude at some distance away 
from the center becomes zero (see the amplitude profile 
at t = t^). After that, the center decays very rapidly (see 
the amplitude profile at t = £5). (Note that in this section 
and the next section, we present numerical results for (0) 
with S r = Si = 0. However, we have also verified that the 
CGL model with parametric forcing also has similar self 
similarity properties.) 



IV. STATISTICS OF BURSTS 

The self similar properties of bursts implies that the 
solutions of the CGL model consist of self-similar bursts 
of various maximum heights. Thus, we expect that the 
enhanced tail (the deviation from the Gaussian distri- 
bution) of the amplitude probability distribution P(|A|) 
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FIG. 8: Statistics of localized events ("bursts"), (a) P hj (\A\) 
at three different times: ti(O) = 10.448 and h = 4.52, 
ts(A) = 89.0 and h = 6.31, and i 3 (D) = 94.0 and h = 3.85. 
(b) The frequency of bursts which have maximum height h, 
g(h). The inset indicates C{h) versus h defined in Eq. (|l3|). 



can be understood by the statistics of bursts. In particu- 
lar, we consider g(h), the frequency of bursts which have 
maximum height h, and a distribution P,-(|^4|) defined for 
an individual burst (burst j), as follows. 

We define the time interval for each burst as the time 
between when its peak value exceeds 2| J 4| a „ g and when its 
peak value drops below 2|A| at)g [typically the time dura- 
tion of a burst is less than 0.2, see Fig. 0(b)]. Here ^4|aug 
is the space average of \A\ over the entire spatial grid of 
the simulation at each time t [\A\ avg is approximately 
constant at about 0.3 over all time steps in the simula- 
tion, see the dashed line in Fig. 0(b)]. Consistent with 
the observation that a typical burst has radial symmetry, 
we define the domain of the burst to be a circular region 
of radius r e // centered at the burst maximum, where 
r e f f is the maximum radius of a circle such that the av- 
erage of \A\ over the perimeter of the circle is greater 
than 2\A\ avg (typically, 1.23 < r eff < 4.91). In Fig. |(a) 
we show the distribution Pj(|A|) for the three bursts in 
Fig. 0(b) (thick solid lines). The first burst (j = 1) 
has h = 4.52 and is plotted as the open circles in Fig. 
||(a); the second burst (j = 2) has h = 6.31 and is plot- 
ted as the open triangles; and the third burst (j = 3) has 
h = 3.85 and is plotted as the open squares. These distri- 
butions are obtained from histograms of the values of \A\ 
at grid points in the domains and time steps in the dura- 
tion of each of these bursts. We obtain g{h) by counting 
the number of bursts which have maximum heights be- 
tween h and h + Ah, where Ah = 0.2 [see Fig. 0(b)]. (In 
Fig. |(a) is not plotted for \A\ < 2|A| n „ fl , since, 

by our procedure this range lacks meaning, and since we 
are interested in the behavior at large values of \A\.) We 
note that the Pj(|.A|) in Fig. ||all approximately coincide 
for \A\ < h. Thus the only characteristic of the bursts on 
which Pj(|j4|) depends is the maximum burst amplitude 
h at which P,-(|A|) goes to zero. To incorporate this fact, 
we replace P,-(|A|) by the notation P^(|A|). 

The above suggests that P(|A|) can be obtained from 
the following approximation [15| 



P{\A\) 



g(h)P h (\A\)dh, 



(11) 
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FIG. 9: P(|A|) versus Circles represent P(|A|) obtained 
directly from our numerical solutions of (^), while pluses rep- 
resent -P(l-Aj) obtained using P*(|A|) and g(h) from Fig. pi 
and Eqs. Q-Q. 

where P/j(|A|) is a probability distribution of a single 
burst whose temporal maximum amplitude is h. Since 
P/j(|A|) vanishes for \A\ > h and because the P^dAQ 
approximately coincide for \A\ < h [see Fig. ||(a)], we 
approximate P^(|^4|) as 

P h (\A\) ~ C-\h)e{h-\A\)P*(\A\), (12) 

C(h) = [ P*QA\)d\A\, (13) 
Jo 

where C(h) is a normalization factor [C(h) ~ 1, when 
h > 1; see the inset on Fig. §(b)], 9(h — \A\) is a step 
function, and P* is the distribution that we numerically 
compute at the largest value of h that we considered 
{hmax — 6.31). Using ( pi] ) and (fl3|), we can further ap- 
proximate P(|A|) as 

poo 

P(\A\) ~ / C-\h)6(h-\A\)g(h)P*QA\)dh 
Jo 

/■oo 
C- 1 (/i).g(/i)d/i. (14) 

(The integral in ( p4| ) is the cumulative frequency of bursts 
which have maximum height greater than \A\.) 

Figures || show the numerically obtained g(h) and 
P*(\A\). Inserting P*(|A|) into Eq.((l|) and Eq.JiJ), we 



obtain the prediction for P(|A|) plotted as pluses in Fig. 
^ for | A | > 2|j4| q , U9 . This appears to agree well with 
the P(|^4|) obtained from our numerical solutions of (||) 
(open circles). (Note that we shift the predicted P(|^4I) 
(pluses) to the P(|A|) (open circles) obtained from (0) 
after removing data points for |^4| < 2|/l| a „ ff .) 



V. CONCLUSION 

We find that the large A behavior of the PDF obtained 
from our CGL solutions is approximately described by 
a stretched exponential form, P(|A|) « e - '" 4 '', where 
i] < 1. In addition, for small A, P r (A r ) and Pi(Ai) are 
approximately Gaussian, as is the case for a random lin- 
ear superposition of waves. We also observe the self simi- 
lar properties of individual bursts, which allow us to con- 
sider the large amplitude behavior of our CGL solutions 
as composites of coherent self similar bursts. Based on 
this we explain the observed non-Gaussian P(|A|) using 
the nonlinear characteristics of individual bursts P/j(|yi|) 
combined with the statistics of burst occurrences g{h). 

These results lead us to conjecture the following pic- 
ture of rare intense events in our model. Linear insta- 
bility leads to a background of relatively low amplitude 
waves that are weakly interacting and result in a random- 
like, incoherent background and low |^4| Gaussian behav- 
ior of P r (A r ) and Pi(Ai). When, by chance, this inco- 
herent behavior results in local conditions conductive to 
burst formation, nonlinear, coherent, self-reinforcing col- 
lapse takes over and promotes a large growth and spiking 
of A. We believe that this general mechanism may be op- 
erative in a variety of physical situations in which rare 
intense events occur {e.g., the water wave experiments of 
Ref. §). 
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events. We thank P. N. Guzdar for his advice on nu- 
merics. This work was supported by the Office of Naval 
Research (Physics) and by the National Science Founda- 
tion (PHY0098632). 
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